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Inspired by Einstein-Podolsky-Rosen-Bohm experiments with photons, we construct an 

C ' 

■ event-based simulation model in which every essential element in the ideal experiment has 



^3 



> 



a counterpart. The model satisfies Einstein's criteria of local causality and does not rely on 
concepts of quantum and probability theory. We consider experiments in which the averages 
correspond to those of a singlet and product state of a system of two S — 1/2 particles. 
The data is analyzed according to the experimental procedure, employing a time window 
^^O ■ to identify pairs. We study how the time window and the passage time of the photons, 

■ which depends on the relative angle between their polarization and the polarizer's direction, 
^SJ ■ influences the correlations, demonstrating that the properties of the optical elements in the 

■ observation stations affect the correlations although the stations are separated spatially and 

■ temporarily. We show that the model can reproduce results which are considered to be 



X 
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intrinsically quantum mechanical. 
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1. Introduction 

Recently, there has been increasing interest in new ways of information processing that 
exploit quantum mechanical correlations. In general, quantum theory describes the state of 
the system by the wavefunction from which we obtain the ensemble averaged value of quanti- 
ties. Quantum theory successfully describes the averaged value of a quantity that is obtained 
in experiments by macroscopic observations. However, when we consider quantum mechanical 
correlations we have to be careful. For example, if we consider the expectation value of a prod- 
uct of two quantities A and i?, we have to measure AB. It is not sufficient to measure A and 
B separately. When we have to measure the correlation AB experimentally, we need a proper 
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definition of the eorrelation. In experiments, we observe aj and fej in tlie i-th. measurement 
for A and B. Likewise, we must properly define the meaning of a "pair" of data {ai,bi) that 
corresponds to AB. If we overlook this point, the interesting nature of quantum correlation 
may disappear. This point is most clearly illustrated by the famous Einstein-Podolsky-Rosen 
(EPR) paradox. 

In 1935, EPR proposed a gedanken experiment, which led them to the conclusion that 
quantum theory is not a complete theory.^ Their reasoning was based on notions about com- 
pleteness, physical reality and locality. Einstein later expressed the principle of locality as 
"T/ie real factual situation of the system S2 is independent of what is done with the system 
Si, which is spatially separated from the former" an ontological definition of locality that 
is now known as Einstein's criteria of local causality. The question arose whether certain 
apparently paradoxical predictions of quantum theory could be experimentally tested. 

Bohm reformulated in 1951 the EPR gedanken experiment into a form which is conceptu- 
ally equivalent but easier to treat mathematically.^ In Bohm's model, a source emits pairs of 
particles with opposite magnetic moments. The two particles separate spatially and propagate 
in free space to an observation station in which they are detected. As the particle arrives at one 
of the two observation stations it passes through a Stern-Gerlach magnet.^ The Stern-Gerlach 
magnet deflects the particle, depending on the orientation of the magnet and the magnetic 
moment of the particle. The deflection defines the spin S = ±1/2 of the particle.*^ As the par- 
ticle leaves the Stern-Gerlach magnet, it generates a signal in one of the two detectors placed 
behind the Stern-Gerlach magnet. The firing of the detector corresponds to a detection event. 

Inspired by Bohm's proposal, Bell derived in 1964 an inequality that imposes restrictions 
on the correlations between the results of the measurements on the two spin-1/2 particles.^ 
Bell demonstrated that the correlation function for the singlet state violates his inequality. 
Hence, quantum theory is in conflict with at least one of the assumptions that were used 
in the derivation of Bell's inequality. The demonstration of the discrepancy between certain 
quantum mechanical expectation values and Bell's inequality is known as Bell's theorem.^ Bell 
concluded that quantum theory is not compatible with Einstein's criteria of local causality 
and that no physical theory of local hidden variables can reproduce all of the predictions of 
quantum theory.^ 

Originally, Bell derived the inequality under the condition that the probability distribution 
of the observation of A is independent of that of B. This is a more restrictive condition than 
the requirement that the physical experimental procedures which axe used to measure A and 
B, are independent. In fact, the inequality may be violated if there is some relation between the 
observations of A and B, regardless whether this relation is of quantum mechanical origin or 
not. In this sense, the popular statement "a classical system cannot violate the Bell inequality" 
is misleading. One has to be very careful and check if the system under study satisfies all the 
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Fig. 1. (color online) Schematic diagram of an EPRB experiment with photons. 

conditions that are necessary to derive the inequahty.^"^ Indeed, it has been pointed out 
that under certain, physicahy reasonable assumptions, a system not relying on any concept 
of quantum theory and that obeys Einstein's criteria of local causality can also violate the 
original Bell inequality.^' ^ The common feature of these models is the presence of a time 
window to identify the single two-particle systems, as in real EPRB experiments. From 
these observations, it is clear that a violation of Bells inequality is not enough to conclude 
that there are quantum correlations.^'^ In this paper, we study the two-particle correlations 
with one of these, what might be called, classical models and we demonstrate that their key 
feature, the dependence of the correlation on the time window that is used to identify the 
pair, allows us to reproduce the correlations that are characteristic for a quantum system of 
two 5 = 1/2 particles. 

For this purpose, we construct a computer model that satisfies Einstein's criteria of local 
causality and use this model to simulate the experiment with optical switches, as performed 
by Aspect et al}^ The sources used in EPRB experiments with photons emit photons with 
opposite but otherwise unpredictable polarization. We refer to this experimental set-up as 
Experiment I. Inserting polarizers between the source and the observation stations changes 
the pair generation procedure such that the two photons have a fixed polarization. We refer 
to this set-up as Experiment II. As a result of the fixed polarization of the photons the photon 
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intensity measured in the detectors behind the polarizers in each observation station obeys 
Mains' law. Our simulation model reproduces the correct quantum mechanical behavior for 
the single-particle and two-particle correlation function for both types of experiments. The 
difference between this model and the model described in Ref.,^ is the algorithm to simulate 
the polarizer. In Ref.^ we used a model for the polarizers that is too simple to correctly 
describe experiments of type II. 

This paper is organized as follows. In Section 2 we describe the experimental set-up, the 
data gathering method and the data analysis procedures used in EPRB experiments with 
photons, and in particular in the timing experiment with optical switches by Aspect et al}^ 
In Appendix A we present an analysis of real experimental data of another EPRB experiment 
with photons^'^ that is the successor of the experiment by Aspect et al}^ We make a distinction 
between experiments of type I and II. A brief review of the analysis of the experiments in the 
framework of quantum theory is given in Section 3. We give explicit expressions for the single- 
and two-particle expectation values for both types of experiments and we introduce Bell's 
inequality. We discuss the fundamental problem of relating quantum theory with the data set 
recorded in the experiment. In Section 4 we describe our computer simulation model of EPRB 
experiments with photons. We give an explicit description of the algorithm to simulate the 
photons, the observation stations containing the polarizers and detectors, and the data analysis 
procedure. Every essential element of the experiment has a counterpart in the algorithm. For 
some model parameters we can compute the two-particle correlation function analytically, as 
shown in Appendix B. In Section 5 wc discuss our simulation and analytical results and where 
appropriate we compare them to the results obtained from quantum theory. Section 6 presents 
a summary and a discussion of our results. 

2. EPRB experiment with photons 

A schematic diagram of the type I timing experiment with optical switches is shown in 
Fig. 1 (see also Fig. 2 in^^). A source emits pairs of photons with opposite polarization. 
Each photon of a pair propagates to an observation station in which it is manipulated and 
detected. The two stations are separated spatially and temporally. This arrangement prevents 
the observation at station 1 (2) to have a causal effect on the data registered at station 2 (1). 

As the photon arrives at station i = 1 and i = 2, it passes through an optical switch 
that directs the photon to one of the two polarizers with a fixed orientation (see Fig. 1). The 
orientation of the two polarizers in each observation station is characterized by the angles 
9i and 9'^. As the photon leaves the polarizer, it generates a signal in one of the two detec- 
tors. Each station has it own clock that assigns a time-tag to each signal generated by one 
of the two detectors. Effectively, this procedure discretizes time in intervals of a width that 
is determined by the time-tag resolution r. The time-tag generators are synchronized before 
each run. This procedure is necessary because in time, the clocks may become unsynchro- 
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nizcd. Furthermore, in real experiments, some photons may not be detected. In this paper, we 
consider ideal experiments only. Hence, our simulation procedure does not allow for practical 
loopholes, such as the detection loophole (lost photons and registration of accidental signals), 
the "fair sampling" loophole, and synchronization problems, that might be present in real 
optical experiments. 

In the experiment, the firing of a detector is regarded as an event. At the nth event, 
the data recorded on a hard disk (not shown) at station i = 1,2 consists of jn,i = ^n.i) ^ 
depending on the state of the optical switches, = ±1, specifying which of the two detectors 
behind the selected polarizer fired and the time tag tn,i indicating the time at which a detector 
fired. Hence, the set of data collected at station i = 1,2 during a run of events may be 
written as 

'^i ~ {Xn,i ~ tn,i-i Tn,i|'^ ~ Ij • • • j • (1) 

Any experimental procedure requires some criterion to decide which detection events are 
to be considered as stemming from a single two-particle system. In EPRB-type experiments 
with photons, this decision is taken on the basis of coincidence in time.^^'^^ Coincidences are 
identified by comparing the time differences {tn,i—tn,2\n = 1, . . . , N} with a time window W }^ 
Thus, for each pair of rotation angles a = Oi^O'-^ and /3 = ^2,^21 the number of coincidences 
between detectors Dx,i {x = ±1), D'^ ^^ {x = ±1) at station 1 and detectors Dy^2 {y = il), 

2 {y = ±1) at station 2 is given by 

N 

Cxy = Cxy{a,P) = Sx,Xn,iSy,Xn,2Sa,'r„,i^/3,'yn.2'^i^ ~ \tn,l - tn,2\), (2) 
n=l 

where &{t) is the Heaviside step function. The correlation E{a, (3) between the coincidence 
counts is then given by 

where the denominator in Eq.(3) is the sum of all coincidences. In practice, the data {Ti, T2} 
are analyzed long after the data has been collected. In general, the numerical values for the 
coincidences Cxy {a, (3) and correlation E{oi, (3) depend on the time-tag resolution and the time 
window used to identify the coincidences. 

In Experiment II, extra polarizers are inserted between the source and the observation 
stations. -"^^ We denote the orientations of these polarizers by the angles r/i and 772- 
2. 1 Role of the time window 

As we already mentioned, in our simulation we leave no room for practical loopholes such as 
detection and "fair sampling" loopholes that may be used to invalidate the conclusions drawn 
from real experiments. The point of view taken in this paper is that we want to perform 
a simulation of ideal experiments and show that we can reproduce the results of quantum 
theory. 
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Most theoretical treatments of the EPRB experiment assume that the correlation, as 
measured in experiment, is given by^ 

N 

^ ~ 5Z ^^'^n,l^y,Xn,2^ (4) 
n=l 

where we assume that the pairs are well defined. This expression, however, is obtained from 
Eq. (2) by taking the limit W — >■ oo, hence the notation ci*^^. An argument that might justify 
taking the limit W ^ oo and hence the expectation that the correlation does not strongly 
depend on W (disregarding statistical fluctuations), is the hypothesis that the time differences 
originate from some random processes that do not depend on the polarization of the photons 
and on the settings of the polarizers. However, the assumption that the time differences are 
independent random variables may not be correct and in fact, in experiments, a lot of effort is 
made to reduce (not increase) W^^ (see also Appendix A). As we will see later, our simulation 
results agree with the results of quantum theory if we assume that the time differences are 
random variables that depend on the settings of the polarizers and the polarization and if we 
consider the limit W ^ 0. 

3. Quantum Theory 

As is well known, quantum theory itself has nothing to say about the individual events as 
they are observed in experiments (quantum measurement paradox), but it provides a frame- 
work to compute the probability for the various possible events to occur. ^' In this section, 
we give a brief account of the quantum mechanical calculation of the averages obtained in 
the EPRB experiment described earlier, strictly staying within the axiomatic framework that 
quantum theory provides. 

In the quantum mechanical description of Experiment I, the source is assumed to emit 
two photons of which the polarization is described by the state 

1^) = -L m,\Vh - \V),\Hh) = {\HV) - \VH)) , (5) 

where H and V denote the horizontal and vertical polarization and the subscripts refer to 
photon 1 and 2, respectively. The state |*) cannot be written as a product of single-photon 
states, hence it is an entangled state. 

In Experiment II, the photons have a definite polarization when they enter the observation 
station. The polarization of the two photons is described by the product state 

I*) =(cos77i|ii")i + sin ?7i I y)i) (cos 7/2 |-H")2 + sin 7721^)2). (6) 

Each of the polarizers in the observation stations splits the beam of incoming photons. 
Using the fact that the two-dimensional vector space with basis vectors {\H), \ V)} is isomor- 
phic to the vector space of spin-1/2 particles, we may use the quantum theory of the latter to 
describe the action of a polarizer as a rotation about its angle of orientation (^i, 6[,62 or 6*2), 
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followed by the measurement of the z-component of the Pauli spin matrix. More specifically, 
a polarizer with orientation a changes the states \H) and \y) according to 

\E.) — cosa|i?) +sina|F), 

|y) ^ -sina|iJ) +cosa|y). (7) 

Hence, the polarizers at station 1 and 2 with orientation a = Q\^Q'y and /3 = respectively, 
change the state |^) into 

\^) = R{a)Rm^), (8) 

where it is implicitly understood that R{ol) and i?(/3) operate on the spin of particle 1 and 2, 
respectively. The rotation matrix -R(0) is given by 

, , / COS0 sin^ \ 

R{e) = _ ) • (9) 

y — sm COS d ) 

According to quantum theory, the expectation value of counting photons at the + (— ) detector 
behind the polarizer with orientation a = Qx^Q'^ {(3 = 92,02) '^s given by^ 

P+{a) = ($|l + c7^|$)/2 = 1/2 + {'^\R-^{a)afR{a)\'^)/2 

= 1/2 + (*|(7iCOs2q; + (7i sin2a|*)/2, 

= (f |1 - al\^)/2 = 1/2 - mR-H(3)alR{P)\^)/2 

= 1/2- (^-laf cos2/3 + c7f sin2/3|*)/2, (10) 

where ai = (erf , cr^, erf ) and (72 = (erf , erf) erf) the Pauli spin- 1/2 matrices for particles 1 
and 2, respectively.^ 

The expectation values of the z-components of the Pauli-spin matrices are given by 

Ei{a) = (^-kf 1$) = P+(q) - P_(q), 

E^iP) = = P+iP) - P^iP). (11) 

The main objective of EPRB experiments is to measure the two-particle correlation 

E{a,P) = ($|(Tfa||$) = {^\R-\a)afR{a)R-\p)alR{p)\^). (12) 

Table I gives the explicit expressions for the expectation values defined by Eqs. (10), (11) 
and (12) for the two different types of experiments. From Table I, it is clear that measuring 
Ei(a), E2{P) and E{a,P) for various a and /3 suffices to distinguish between systems in the 
entangled state (Experiment I) or in the product state (Experiment II). 

Data of EPRB experiments are often analyzed in terms of the function^^' 

S{a, a', P, P') =E{a, P) - E{a, p') + £;(a', P) + E{a' , /?'), (13) 
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Table I. The single- and two-particle expectation values defined by Eqs. (10), (11) and (12) for the 
two experiments described by the states Eqs. (5) and (6), respectively. 

Experiment I Experiment II 
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m) 
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cos2(a — 
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cos2(/?- 


m) 


E{a,l3) 


— cos 2(a — (3) 


cos 2{a — r)i) cos 2(/? — 


m) 



because it provides clear evidence that the system is described by an entangled state. The 
idea behind this reasoning is that for any product state 

-2 < 5(a,a',/3,/3') < 2, (14) 

an inequality known as one of Bell's generalized inequalities.^*^ This can be seen as follows. For 
any product state |^'), we have E{a,(3) = Ei{a)E2{(3). Let us denote a = Ei{a), b = Ei(a'), 
c = E2{P), and d = E2{[3'). Clearly, a, 6, c, d G [—1, 1]. For any a, 6, c, d G [—1, 1] we have^"^ 

\ac — ad + bc + bd\ < \ac — ad\ + |6c + M| < \a\ \c — d\ + \b\\c + d\ 

<\c-d\ + \c + d\<l-cd+l + cd 

< 2, (15) 

hence Eq. (14) follows. Thus, we conclude that if 1^*) can be written as a product state, we 
must have 

Smax= max \S{a,a',/3,l3')\<2. (16) 

a,a',l3,l3' 

Furthermore, it can be shown that^^ 

\S{a,a',l3,l3')\<2V2, (17) 

independent of the choice of l^*). In other words, if 2 < Smax < 2\/2, the quantum system is 
in an entangled state. 

For later use, we introduce the function 

S{e) = S{a,a + 2e,a + e,a + 3e), (18) 

where we have fixed the relation between the angles /? = a+9, a' = a +26, f3' = a+39 through 
the angle 6. Because of rotational invariance, S{6) docs not depend on a and therefore, we 
set a = to simplify matters a little. In the case of Experiment I, E{a,P) = — cos2(a — P) 
and we find 

5(0) = 3 cos 20 -cos 6^, (19) 
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which reaches its S{9) = 2V2 at 9 = Tr/8 + j7r/2, where J IS an 

integer number. 

Analysis of the experimental data/'^"^'^' yields results that are in good agreement with 
the expressions in Table I, leading to the conclusion that in a quantum mechanical description 
of Experiment I, the state does not factorize, in spite of the fact that the particles are spatially 
and temporally separated and do not interact. Our analysis of the few experimental data for 
Experiment I that is publicly available^^ supports this conclusion (see Appendix A). 
3.1 From quantum theory to data 

According to the formalism of quantum theory, the result of a measurement is an eigen- 
value of the dynamical variable that is being measured.^ Applied to the case of the EPRB 
experiment, each measurement yields an eigenvalue of the matrices A = R~'^{a)afR{a), 
B = R-'^{(3)a^R{(3), and C = R-\a)R-^(3)ala^R{f3)R{a). Obviously, the nth measure- 
ment of A, B, or C yields an eigenvalue a„ = ±1, = ±1, or c„ = ±1, respectively. 

The conventional interpretation of quantum theory asserts that the outcome of each mea- 
surement constitutes a Bernoulli trial, that is we assign the same probability to an outcome, 
independent of which trial is considered and independent of what happened in any of the 
other measurements. In other words, the probability to observe for instance is logically 
independent from the probability to observe for all n ^ m. 

For simplicity, we now focus on the case where all photons arc directed towards the po- 
larizers with orientation 9i. Let us then inquire how we can simulate the quantum mechanical 
results of the EPRB experiment (sec Table I) without leaving the framework of quantum 
theory. Evidently, this is a nearly trivial exercise. All we have to do is set up three Bernoulli 
processes that generate sets of data {an = ±1, bn = ±1, Cn = ±l|n = 1, . . . , N} such that 
TV ^ N ^ N 

-J^^'ri^E.iei), -J2bn^E2{e2), - C„ « ^(^1 , ^2) , (20) 

n=l n=l n=l 

for all 9i and large N. However, this line of reasoning brings out the fundamental problem of 
relating the set of data 

Q = {an,bn,Cn\n = l,...,N}, (21) 
obtained from quantum theory with the set of data 

£ = {Xn,l,Xn,2, tn,l,tn,2\l^ = !,•••, N}, (22) 

recorded in the experiment. 

In the case that we measure a property of a single particle, we may identify Xn,i with a„ 
and Xn,2 with bn, yielding 

N N 

Ei{9i), ~ E2{92), (23) 

n=l n=l 
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where the prime indicates that the sum runs over all events that yield a coincidence, that is 
for all events for which Q{W — \tn,i — tn,2\) = 1- 

As we know from the work of Bell and others,^' simple assignments of the form c^i ^— 
Xn,iXn,2 Cannot reproduce the result of quantum theory for E{9i,62)- However, as we need 
to identify pairs for measurements that involve properties of two particles, there is no simple 
a-priori rule to relate c„ to the data {xn,i, Xn,2i tn,i,'tn,2}- In EPRB experiments with photons 
the correlation £'(^1,^2) is calculated according to Eq. (3), using the coincidences Eq. (2). 
That is, we adopt the assignment 

^ -n,lXr,2e{W-K^-t^,2\)^^ (24) 

En=l®iW-Kl-tn,2\) 

In Eq. (24), the coincidence window W enters because it is necessary to have a criterion to 
decide which particles belong to a single two-particle system, an essential ingredient in any 
real EPRB experiment. As this choice is, in a sense, ad hoc, we have to study the significancy 
and implications of this choice, as we will do in the next sections. 

4. Simulation model 

We now take up the main challenge, the construction of processes that generate the data 
sets Eq. (1) such that they reproduce the results of quantum theory, summarized in Table I. A 
concrete simulation model of the EPRB experiment sketched in Fig. 1 requires a specification 
of the information carried by the particles, of the algorithm that simulates the source and the 
observation stations, and of the procedure to analyze the data. From the specification of the 
algorithm, it will be clear that it complies with Einstein's criteria of local causality on the 
ontological level: Once the particles leave the source, an action at observation station 1 (2) 
can, in no way, have a causal effect on the outcome of the measurement at observation station 
2(1). 

4-1 Source and particles 

The source emits particles that carry a vector Sn,i = (cos(^„ + ~ l)7i"/2), sin(^„ + ~ 
l)7r/2)), representing the polarization of the photons. The "polarization state" of a particle is 
completely characterized by which is distributed uniformly over the interval [0, 27r[. We use 
uniform random numbers to mimic the apparent unpredictability of the experimental data. 
However, from the description of the algorithm, it trivially follows that instead of uniform 
random number generators, simple counters that sample the interval [0, 27r[ in a systematic, 
but uniform, manner might be employed as well. This is akin to performing integrals by the 
trapezium rule instead of by Monte Carlo sampling. 
4-2 Observation stations 

The input-output relation of a polarizer is rather simple: For each input event, the al- 
gorithm maps the input vector S onto a single output bit x. The value of the output bit 
depends on the orientation of the polarizer a = (cos a, sin a). According to Malus' law, for 
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fixed S = (cos ^, sin ^) and fixed a, tfie bits to be generated sucli tliat 



— ^Xn = cos2{^-a), 



(25) 



n=l 



with probability one. If, as in Experiment I, the input vectors S are distributed uniformly 
over the unit circle, the sequence of output bits should satisfy 



with probability one, independent of the orientation a of the polarizer. 

As we work under the hypothesis of ideal experiments, the algorithm to simulate each of 
the four different polarizers in Fig. 1 should be identical. Evidently, for the present purpose, 
if we switch from Experiment I to Experiment II, it is not permitted to change the algorithm 
for the polarizer. 

In this paper, we use a deterministic model for a polarizer. Elsewhere, we have demon- 
strated that simple deterministic, local, causal and classical processes that have a primi- 
tive form of learning capability can be used to simulate quantum systems, not by solving a 
wave equation but directly through event-by-event simulation.^^^'^ The events are generated 
such that their frequencies of occurrence agree with the quantum mechanical probabilities. In 
this simulation approach, the basic processing unit is called a deterministic learning machine 
(DLM).^^^'^ A DLM learns by processing successive events but does not store the data con- 
tained in the individual events. Connecting the input of a DLM to the output of another DLM 
yields a locally connected network of DLMs. A DLM within the network locally processes the 
data contained in an event and responds by sending a message that may be used as input for 
another DLM. DLMs process messages in a sequential manner and only communicate with 
each other by message passing. In a simple physical picture, a DLM is a device that exchanges 
information with the particles that pass through it. It learns by comparing the message car- 
ried by an event with predictions based on the knowledge acquired by the DLM during the 
processing of previous events. The DLM tries to do this in an efficient manner, effectively by 
minimizing the difference of the data in the message and the DLM's internal representation 
of it?^~^^ We now describe the DLM that simulates the operation of a polarizer. 

Let us focus on the polarizer with orientation 9i . The DLM has an internal two-dimensional 
unit vector R„ = (x„, y„) and a parameter < Z < 1 that controls the pace of learning (to be 
discussed later). The orientation a = (cos ^i, sin ^i) of the polarizer is also part of the input 
to the DLM. The DLM receives as input, the sequence of unit vectors S„ = (cos sin ^„) 
for n = 1, . . . , A?^ where N is the total number of events, generated by the source. We know 
that the output signal of a polarizer depends on the difference between the polarization and 
the orientation of the polarizer only (Mains' law). This is taken into account by rotating S„ 
about 6i. We denote the resulting vector by Y„ = (cos(^„ — ^i), sin(^„ — ^i)). For each input 




(26) 



n=l 
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event, the DLM computes eight trial vectors according to the following rules^^'^° 

Xji = Is Xji 2 1" S -yl — Z^+Z^E^ 2 ) 

y„ = Is'yn^^ + s^l-P + l^yl^^, (27) 

where s,,s' = ±1 are variables that allow us to generate the trial vectors in each of the four 
quadrants and A = ±1 determines whether the trial vector is obtained by rcscaling the x- 
coordinate (A = 1) or the y-coordinate (A = — 1). Note that each of the eight rules generates 
a unit trial vector. 

The final step in the DLM algorithm consists of comparing the vector Yn with each of 
the eight trial vectors {xn,yn)- The DLM updates its internal vector R„ by choosing the trial 
vector, that is the triple (A, s, s'), for which the Euclidian distance || Y„ — y„) || is minimal. 
If the trial vector was obtained by applying a A = +1 rule, the DLM generates a +1 output 
event. Otherwise it generates a —1 output event. 

Elsewhere, we have shown by means of simulations and analytical methods that the DLM 
described above generates —1 and +1 events, that are distributed according to Malus' law 
if the input vector S„ = (cos sin does not change during a sufficiently long sequence 
of events. ^^'^'^ If the input vector S„ = (cos^,^, sin^„) is uniformly distributed over the unit 
circle, the dynamics of the DLM generates events according to the function sign(cos2(^„ — ^i)) 
(results not shown). 

The DLM is a machine with elementary learning capabilities: It is an adaptive system that 
learns from the input events.^^'^^'^*^ The parameter < / < 1 controls the speed of the learning 
process and the accuracy with which the internal vector can represent input vectors. If / is 
close to one, the DLM learns slow and gives an accurate response but it also "forgets" slow, 
that is if the DLM is offered different input vectors, it may take long before its internal vector 
has adapted to the new situation. In a sense, I controls the "coherence" of the system.^^'^'''^'' 
As the DLM, operating according to the rules Eq. (27), can reproduce both Mains' law and 
the function sign(cos2(^„ — Oi)), all that is left to do to completely specify the model of the 
polarizer is to add the mechanics for the time tagging. 

To assign a time-tag to each event, we assume that as a particle passes through the 
detection system, it may experience a time delay. This is a key assumption in the construction 
of the simulation model. In principle, the time-delay of the individual photons cannot be 
derived from the Maxwell equations because they describe waves, not particles. Thus, to 
find an a-priori justification for the assumption that the particle experiences a time-delay 
we are limited to making inferences from experimental data. Empirical evidence is provided 
in Appendix A where we analyze experimental data of an EPRB experiment with photons 
and demonstrate that the average time-of-flight of the photons depends on the orientation 
of the polarizer. Thus, there is experimental evidence that supports the assumption that the 
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photons experience a time delay as they pass through the polarizers. It is unfortunate that the 
experimental data that is publicly available is far too scarce to allow a detailed analysis of the 
time-delay mechanism. Therefore, as a model of the time delay, we will choose a specific model 
that is as simple as possible, is in concert with empirical knowledge, is capable of reproducing 
the results of quantum theory and allows an analytical treatment in particular limiting cases. 

In our model of the time delay, tn,i for a particle is assumed to be distributed uniformly 
over the interval [to,tQ +T]. In practice, we use uniform random numbers to generate tn.i- As 
in the case of the angles the random choice of tn,i is merely convenient, not essential. From 
Eq.(2), it follows that only differences of time delays matter. Hence, we may put to = 0. The 
time-tag for the event n is then G [0, T]. 

We now come to the point that we have to specify T explicitly. In fact, there are not many 
options. Is is an experimental fact that the output intensities of the two light beams emerging 
from a polarizer depend on the direction of polarization of the incident light, relative to the 
angle of the main optical axis of the polarizer crystal. "^-'^ Thus, at least macroscopically, the 
system consisting of a polarized light beam and polarizer is invariant for rotations about the 
direction of propagation of the light wave. Assuming that this invariance carries over to the 
individual particles, we can construct only one number that depends on the relative angle: 
Sn-a, implying that T = T(^„ — can depend on — only. Furthermore, consistency with 
classical electrodynamics requires that functions that depend on the polarization have period 
vr.^^ Thus, we must have T(^„ — 6i) = F((S„^i • a)^). Of course, the arguments that have 
been used to arrive at this form are not sufficient to uniquely fix the form of _F((S„^i • a)^). 
As explained earlier, the available experimental data does not suffice to determine the form 
of F((S,„.i ■ a)'^). Therefore we use simplicity as a criterion to select a specific form. By trial 
and error, we found that T(^„ - ^i) = ToF(|sin2(^„ - 0i)|) = ro|sin2(^„ - Oi)]'^ yields 
useful results. Here, Tq = maxg T{6) is the maximum time delay and defines the unit of time, 
used in the simulation and d is a free parameter of the model. In our numerical work, we 
set To = 1. As we demonstrate later, our model reproduces the quantum results of Table 
I under the hypothesis that the time tags tn,i are distributed uniformly over the interval 
[0, 1 sin2(^„ — ^i)!*^] with d = 2. Needless to say, we do not claim that our choice is the only 
one that reproduces the results of quantum theory for the EPRB experiments. 

A last point left to address on the algorithm for the polarizer is the deterministic character 
of the output sequences generated by the DLMs. In fact, it is easy to change this deterministic 
process into a random process without changing the distribution of and —1 events. A simple 
way to do this is to attach to the particles, a random bit 2„ = ±1, generated by the source. 
Then, we modify the algorithm for the polarizer such that instead of selecting output channel 
A, it selects output channel 2;„A. This procedure has no effect on the DLM dynamics and 
the single- and two-particle counts but, it renders the outputs of the detectors unpredictable. 
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Fig. 2. (color online) Left: Computer simulation of Experiment I in which the source emits photons 
with opposite random polarization (EPRB experiment). Right: Computer simulation of Experi- 
ment II in which the source emits photons with fixed polarization. Solid circles (red): Simulation 
results using DLMs with the time-delay mechanism {d = 2) for the polarizers. Open circles (black): 
Simulation results using DLMs but without using the time-tags (equivalent to d = or cx)) 
to compute the two-particle correlation. Other markers: Average single-particle counts on the 
detectors (see Fig. 1). Squares (green): ^+(6*1) = P+{9[)] Diamonds (green): P-{9i) = P^{9[); 
Plusses (blue): ^+(6*2) P+{0'2); Crosses (blue): ^-(6(2) = P-{0'2). In Experiment I (left), these 
four symbols lie on top of each other. In Experiment II (right), these markers show the typical 
Malus law behavior. Solid line: Quantum theory for SiO). Dashed line at S{9) = 2\/2: Maximum 
of S{9) if the system is described by quantum theory. Dashed line at S{9) = 2: Maximum of S{9) 
if the system is described by the class of models introduced by Bell;^ Dashed hue at S[9) = 1/2: 
Expected number of -1-1 and —1 events recorded by the detectors if the input to the polarizers 
consist of photons with random polarization. Dotted lines: Quantum theory for P^{9i) = Pj^{9[), 
P.{9,) = P_{9[), P+{92) = P+(0i) and P_(02) = P-(^^). 

This technical finesse does not affect the final result for the expectation values and therefore 
we disregard it in what follows. 
4.3 Data analysis 

For fixed N and fixed angles 6i and 0- of the polarizers in the observation stations, the 
algorithm described earher generates the data sets Tj, just as experiment does. In order to 
count the coincidences, we choose a time-tag resolution r and a time window W such that 
< T < Tq and r < W. We set the single particle counts Px{a), Py{(3) and the coincidence 
counts Cxy{oi,(3) with a = 0i,9[ and (3 = 02,02: to zero for all x,y = ±1. Then, we make 
a loop over all events Xn,i = ±1 in the data sets and we read off ^n,i = 0n,i,0'ni- To count 
the coincidences, we first compute the discretized time tags kn,i = \tn,i/T'\ for all events in 
both data sets. Here \x\ denotes the smallest integer that is larger or equal to x, that is 
[x] — 1 < X < \x\. According to the procedure adopted in the experiment, an entangled 
photon pair is observed if and only if |/c„_i — kn,2\ < k = [VF/r]. Thus, if \kn^i — kn^2\ < k, we 
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Fig. 3. (color online) Left: Same as Fig. 2 (left) except that d = 1. Dotted line: S{9) calculated ft-om 
Eq.(33). Right: Same as Fig. 2 (left) except that d = 4. Dotted line: S{6) calculated from Eq.(34). 

increment the count Cx„,i,x„,2(cK) /?) and we increment the corresponding single particle counts 

We emphasize that the simulation procedure counts all events that, according to the same 
criterion as the one employed in experiment, correspond to the detection of single two-particle 
systems. 

5. Simulation results 

We use the computer model, described earlier to simulate the experiment depicted in 
Fig. 1. Each polarizer in Fig. 1 is simulated by the same algorithm. The procedure to direct 
the particles to the polarizers is the same as in the laboratory experiment. For each particle 
that enters station 1 (2), a random number generator at station 1 (2) determines which of 
the two polarizers will receive the particle. In practice, we use two different random number 
generators for station 1 and 2 (we have never seen any statistically significant effect of using 
the same one for both stations). The source always sends out two particles with orthogonal, 
random polarization. Unless we insert additional polarizers between the source and stations 1 
and 2, this setup simulates Experiment I. Inserting additional polarizers between the source 
and stations 1 and 2 is the same as sending the particles with fixed polarization to stations 1 
and 2. In this case, we simulate Experiment II. 

The simulation proceeds in the same way as in the experiment, that is we first collect the 
data sets Ti and T2 for various settings of the polarizers (various 7n.i), and then compute 
the single particle counts Eq. (10), the coincidences Eq.(2) and the correlation Eq.(3), from 
which we can calculate the function S{6) (see Eq.(18)). 

In Fig. 2 (left), we present our simulation data for Experiment I, that is for the case that 
the source emits particles with an opposite, random polarization, corresponding to the singlet 
state in the quantum mechanical description. The parameters in these simulations are k = 1, 
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d = 0,2, T = 0.00025, / = 0.999, and N = 10^. The results are not sensitive to the choice of 
these parameters. For instance, the figures (not shown) with the results for k = 1, d = 0,2, 
T = 0.25, / = 0.999, and N = 10^ are barely distinguishable from those of Fig. 2 and, as 
can be expected on general grounds, increasing the number of events N simply reduces the 
fluctuations. The parameter I controls the accuracy with which we can resolve differences in 
the angles: The closer / is to one, the higher the accuracy. We have chosen I such that, with 
the resolution used to plot the data, the effect of / on the results cannot be noticed. The main 
reason for showing the data for = 1, d = 0, 2, r = 0.00025, / = 0.999, and N = 10^ is that 
it allows us to demonstrate that our simulation results are in excellent agreement with the 
analytical results of Appendix B. 

It is clear that for d = 2, the simulation model reproduces the results of quantum theory for 
the single-particle expectation values P±{a) and P±{(3) (see Table I) and S{9) (see Eq. (19)). 
Indeed, the frequency with which each detector fires is approximately one-half and S{9) agrees 
with the expression Eq. (19) for the singlet state. Also shown in Fig. 2 (left) are the results for 
S{9) if we disable the time-delay mechanism. Effectively, this is the same as letting the time 
window ^ oo or setting d = 0. Then, our simulation model generates data that satisfies 
|>S'(^)| < 2, which is what we expect for the class of models studied by Bell.^ 

In Experiment II, the source emits particles with a fixed (not necessarily opposite) polar- 
ization. In the right panel of Fig. 2, we present results, obtained by the same simulation algo- 
rithm as for Experiment I, for the case 9i = O't^ = a = a' = 9 and 92 = 92 = fi = (3' = 9 + 7r/4. 
The angle ^ of the particles is 7r/6 (corresponding to r/i = 7r/6 and r/2 = '7r/6-|-7r/2 in 
the quantum mechanical description). For this choice, we have P+{a) = cos^(0 — "^/G), 
P+{f5) = cos2(^-7r/6-7r/4), E{a,f5) = 2"! sin4(7r/6 - ^) and S{9) = sin4(7r/6 - 0). Except 
for the properties of the particles, the model parameters for Experiment I and II are the same. 

From Fig. 2, it is clear that the event-by-event simulation reproduces the single- and two- 
particle results of quantum theory for both Experiment I and II, without any change to the 
algorithm that simulates the polarizers. 

Having established that the data generated by our "non-quantum" system agrees with 
quantum theory, it is of interest to explore if these dynamical, adaptive systems can generate 
data that is not described by quantum theory or by the simple, locally causal probabilistic 
models introduced by Bell.^ We can readily give an affirmative answer to this question by 
repeating the simulations for Experiment I (see Fig. 2 (left)) for different values of the time- 
delay parameter d, all other parameters being the same as those used to obtain the data 
presented in Fig. 2. 

For d = 0, simulations with or without time-delay mechanism yields data that, within 
the usual statistical errors, are the same (results not shown) and satisfy 15(^)1 < 2 . Figure 3 
shows the simulation data for d = 1 and d = 4. For < d < 2 our model yields two-particle 
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Fig. 4. (color online) Smax = ^8.Xa,a' .p,i3' S{a, a' , j3, j3') as a function of the time window W relative 
to the time-tag resolution t. Curves from bottom to top: Results for d = 0, 1, . . . , 10. Dashed line: 
Value of Smax — 2\/2 if the system is described by quantum theory. 

correlations that are stronger than those of the Bell-type models but they are weaker than in 
the case of the singlet state in quantum theory. Therefore, the maximum of S{9) is less than 
2^/2 but larger than two. For c? > 3, we find that the two-particle correlations are significantly 
stronger than in the case of the singlet state in quantum theory. From Fig. 3 it can be seen 
that for d = \ and d = 4 there is good agreement between the results obtained with our 
event-based simulation model and the analytical result for |5'(^)| obtained from Eqs. (33) and 
(34), respectively (see Appendix B). For d = 1, the simulation results show larger fluctuations 
than for d = 4, but in all cases they can be reduced by increasing N (results not shown). 

The simulation results presented in Figs. 2 and 3 have been obtained for W/t = 1 and 
small T (recall that the unit of time in our numerical work is set equal to one). In general, 
in experiment the two-particle correlation depends on both W and r. Our simulation model 
makes definite predictions for this dependence. This can be seen from Fig. 4 which shows 
Smax = ^SiXa^a' ,i3,f3' S {a, a' , (3 , (3') as a function of W/t for various values of d. Smax is cal- 
culated numerically using Eqs. (28) and (29) (see Appendix B). The numerical results agree 
with the values of Smax that have been obtained analytically for W = t ^ 0, d = 0,2 and 
W CO. For d < 2, 2 < Smax < 2\/2 for any value of W/t. Hence, for d < 2 our model 
cannot produce the correlations of the singlet state. For d = 2, 2 < Smax < 2-v/2 and our 
model produces the correlations of the singlet state if W/t — > 0. For d > 2, 2 < Smax < 4, and 
for a range of W/t, Smax > 2\f2, implying that our model exhibits correlations that cannot 
be described by the quantum theory of two spin-1/2 particles. 
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6. Discussion 

Wc have pointed out that a violation of Bell's inequality is not an absolute criterion for 
having quantum correlations. The fact that our event-based model can produce correlations 
that violate the Bell inequality, by itself, is not a surprise, because the expression for the co- 
incidences Eq. (2) cannot be written in a form that allows a derivation of the Bell inequality.^ 
Given the two data sets, recorded at the two observation stations as described in Section 2, 
we may find correlations between the data of both sets that may or may not violate the Bell 
inequality: Mathematically such correlations cannot be excluded.^ Our results do not contra- 
dict the folklore about Bell's theorem. Bell's notion of locality is an attempt to incorporate 
Einstein's criteria of local causality on the ontological level in probabilistic theories. In fact, 
in the derivation of Bell's inequality a strong mathematical (in terms of probability theory) 
assumption was made, namely that any logical relation between both data sets is prohibited. 
Physically, this assumption is usually erronously associated with the independence of opera- 
tions at distant positions. Thus the presence of any physical relation, which could be classical 
or quantum in origin, can lead to a violation of the bounds in Bell's inequality. 

The time-delay mechanism is an example of such a physical relation. In our event-based 
model, the expression for the coincidences is the key ingredient to reproduce the quantum 
mechanical results for the two-particle correlation of the EPRB experiment. This expression, 
based on the time tags of the detection events, is the same as the one employed in EPRB 
laboratory experiments with photons. With this example, we may say that a simulation model 
that strictly satisfies Einstein's criteria of local causality can reproduce the quantum theo- 
retical results for EPRB experiments, without using any concept from quantum theory. That 
is, although our event generating and measurement processes are of classical origin, they still 
lead to a violation of the original Bell inequality. This fact must be considered seriously in any 
real experiment. Instead of the time-delay, we may consider various other mechanisms which 
can cause a similar effect. In fact, there may be many "unpaired signals". This means that 
many photons are destroyed by some reason. So far the destruction is considered to happen 
randomly. But if it would depend on the relative angle between the photon polarization an the 
direction of the polarizer, similar effects as those studied in the present paper could occur. In 
any case, we have to be careful in eliminating those possibilities while studying the appearance 
of quantum correlations. The dependence on the time window is a good check of this fact. 

We have shown that in our model, in the case of Experiment I, the two-spin correlation 
depends on the value of the time window W. By reducing W from infinity to zero, this 
correlation changes from typical Bell-like to singlet-like, without changing the procedure by 
which the particles are emitted by the source. Thus, the character of the correlation not only 
depends on the whole experimental setup but also on the way the data analysis is carried out. 
Hence, from the two-spin correlation itself, one cannot make any definite statement about the 
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character of the source. Thus, the spin correlation is a property of the whole system (which 
is what quantum theory describes), not a property of the source itself. It is of interest to note 
that if we perform a simulation of Experiment II the single-spin and two-spin correlations 
do not depend on the value of the time window W. In this case, the observation stations 
always receive particles with the same polarization and although the number of coincidences 
decreases with W (and the statistical errors increase), the functional form of the correlation 
does not depend on W. 

We have also presented a rigorous proof that our simulation model reproduces the two-spin 
correlation that is characteristic for the singlet state. Furthermore, our model also allows us 
to explore phenomena that cannot be described by quantum theory of two S = 1/2 particles. 

Finally, we also examined the effect that the pair identification criterion has on the two- 
particle correlations for a set of experimental data that is publicly available. The results, 
presented in Appendix A, show a tendency that is similar to the predictions of our simulation 
model, namely that the time window, used as a criterion to identify photon pairs based on 
the time-tag data of single photon events, should be chosen as small as possible in order to 
find results for the single-particle counts and two-particle correlations that agree with the 
quantum theoretical expectation values for a system of two S = 1/2 particles. 

Appendix A 

We illustrate the importance of the choice of the time window W by analyzing a data set 
(the archives Alice.zip and Bob.zip) of an EPRB experiment with photons that is publicly 
available.^^ Technically, the experiment of Ref.^^'^^ is different from the one sketched in Fig. 1, 
but conceptually both experiments are the same. The data in the archives Alice.zip and 
Bob.zip are data for Experiment I in which 9i = 0, 9[ = 7r/4, 62 = 7r/8, and 62 = Stt/S. 

In the real experiment, the number of events detected at station 1 is unlikely to be the same 
as the number of events detected at station 2. In fact, the data sets of Ref.^^ show that station 
1 (Alice.zip) recorded Ni = 388455 events while station 2 (Bob.zip) recorded N2 = 302271 
events. Furthermore, in the real EPRB experiment, there may be an unknown shift A (assumed 
to be constant during the experiment) between the times {tn,i\n = 1, • • • ,iVi} gathered at 
station 1 and the times {tm,2\'m = 1, • • • ,N2} recorded at station 2. Therefore, there is some 
extra ambiguity in matching the data of station 1 to the data of station 2. 

A simple data processing procedure that resolves this ambiguity consists of two steps.^^ 
First, we make a histogram of the time differences tn,i — im,2 with a small but reasonable 
resolution (we used 0.5 ns). Then, we fix the value of the time-shift A by searching for the 
time difference for which the histogram reaches its maximum, that is we maximize the number 
of coincidences by a suitable choice of A. For the case at hand, we find A = 4 ns. Finally, we 
compute the coincidences, the two-particle average, and Smax using the same expressions as the 
ones used to analyze the computer simulation data. The average times between two detection 
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Fig. 5. (color online) Smax as a function of the time window W, computed from the data sets con- 
tained in the archives Alice.zip and Bob. zip that can be downloaded from Ref.^^ Bullets (red): 
Data obtained by using the relative time shift A = 4 ns that maximizes the number of coinci- 
dences. The maximum value of Smax ~ 2.73 is found sX W — 2 ns. Crosses (blue): Raw data 
(A = 0). The maximum value of Smax ~ 2.89 is found at = 3 ns. Dashed line at \S{9)\ — 2a/2: 
Smax if the system is described by quantum theory (see Section 3). Dashed line at \S{9)\ — 2: 
Smax if the system is described by the class of models introduced by Bell.^ 

events is 2.5 ms and 3.3 ms for Alice and Bob, respectively. The number of coincidences (with 
double counts removed) is 13975 and 2899 for (A = 4 ns, = 2 ns) and {A = , W = 3 ns) 
respectively. 

In Fig. 5 we present the results for Smax as a function of the time window W. First, it is 
clear that Smax decreases as W increases. Second, the procedure of maximizing the coincidence 
count by varying A reduces the maximum value of Smax from a value 2.89 (A = 0) that 
considerably exceeds the maximum for the quantum system (2-v/2, see Section 3) to a value 
2.73 (the value cited in Ref.^^) that violates the Bell inequality and is less than the maximum 
for the quantum system. The optimized experimental results (bullets in Fig. 5) and the results 
of our simulation model (see Fig. 4, third line from the bottom ) are qualitatively very similar. 

The fact that the "uncorrected" data (A = 0) violate the rigorous bound for the quantum 
system should not been taken as evidence that quantum theory is "wrong" : It merely indicates 
that the way in which the data of the two stations has been grouped in two-particle events is 
not optimal. 

Finally, we use the experimental data to show that the time delays depend on the orienta- 
tion of the polarizer. To this end, we select all coincidences between -D+,i and (see Fig. 1) 
and make a histogram of the coincidence counts as a function of the time-tag difference, for 
fixed orientation = and the two orientations O2 = vr/8, 37r/8 (other combinations give sim- 
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Fig. 6. (color online) Normalized coincidence counts as a function of time tag difference i„_i — tn^2, 
computed from the data sets contained in the archives Alice.zip and Bob. zip, using the relative 
time shift A = 4 ns that maximizes the number of coincidences. Bullets (red): 9i — and 62 — 7r/8; 
Crosses (blue): di — and 62 — Stt/S. 

ilar results). The results of this analysis are shown in Fig. 6. The maximum of the distribution 
shifts by approximately 1 ns as the polarizer at station 2 is rotated by -7r/4, a demonstration 
that the time-tag data is sensitive to the orientation of the polarizer at station 2. A similar 
distribution of time-delays (of about the same width) was also observed in a much older ex- 
perimental realization of the EPRB experiment. The birefringent properties of the optical 
elements (polarizers and electro-optic modulators) might be responsible for this time delay. 
A more detailed quantitative and exploratory analysis of this time delay requires dedicated 
retardation measurements for these specific optical elements in single-photon set-ups. 

Appendix B 

In the case of Experiment I and for some choices of the model parameters, we can compute 
the correlation Eq. (3) analytically.^ In the limit N ^ 00, Eq. (3) can be written as 



where P{Ti,T2,W) is the density of coincidences for fixed (a, (3) and polarization angle 
(within a small interval d^), Ti = \ sin 2(^ — a)!*^, T2 = \ sin 2(^-/5) j*^, xi(^,a) = sign(cos2(^ — 
a)), and X2{^, 13) = sign(cos 2(^ - /5)). 

Mathematically, expression Eq. (28) can give rise to almost any correlation. The density 
of coincidences appears as a result of integrating over the distribution of time tags and it is 
this integration that leads to the appearance of the correlations. Of course, the specific form 
of the correlations depends on the choice of the functional dependence of Ti and T2 on the 




(28) 
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Fig. 7. (color online) Graphical representation of the process of counting pairs. The time interval 
is divided in bins of size r, represented by the elementary squares. The two parallel, 45° lines 
indicate the time window W, which was chosen to be 2t in this example. In the limit N oo, the 
total number of pairs for fixed a, f3, and ^ is given by the number of whole squares that fall within 
the time window and satisfy 1 < ki < Ki for i = 1,2. For Ki > K2, all filled squares contribute 
while for K[ — K2, the dark gray square does not contribute. For Ki < K2 we interchange labels 
1 and 2. 



angles. 

The expression for P{Ti,T2, W) can be derived as follows. For a fixed time-tag resolution 
< r < 1, the discretized time-tag for the nth detection event is given by kn^i = \tn,iT~^~\ 
where [x] denotes the smallest integer that is larger or equal to x. The discretized time-tag 
kn,i takes integer values between 1 and Ki = [r~^Tj] , where Ki is the maximum, discretized 
time delay for a particle with polarization ^ and passing through the polarizer with orientation 
7j, where 71 = a and 72 = (3. If \kn^i — kn,2\ < k = \t~^W~\ , the two photons are defined to 
form a pair. For fixed a, /?, and ^, we can count the total number of pairs, or coincidences 
C, by considering the graphical representation shown in Fig. 7. After a careful examination 
of all possibilities, we find that 

C = C{Ki, K2, k) = {2ko - l)ki2 - koiko - l)/2 

- max(0, {K12 - 1) max(0, i^i2)/2) 

-|- max(0, k — ko)ko 
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- max(0, kku - K1K2), (29) 

where ko = mm{Ki, K2, k), k\2 = inin(ifi, ^2), and K\2 = ki2 — max(0, m8ix{Ki, K2) — k). 

It is clear that the result for the coincidences depends on the time-tag resolution r, the 

time window W and the number of events N, just as in real experiments. ^'^"^^'^^ Formula 

Eq. (29) greatly simplifies if we consider the case k = 1 {W = r), yielding C{Ki,K2,l) = 

m.m{Ki, K2) as is evident by looking at Fig. 7. For fixed a, /3, and ^, and W = t, the density 

P{Ti,T2,t) = C{Ki,K2,l)/KiK2 that we register two particles with a time-tag difference 

less than r is bounded by 

minfTi + r, T2 + r) „ , min(Ti,T2) , , 

^ (T LvT 1^/ < P{Ti,T2,t) < T (30) 

For = T ^ and d = 2, Eq. (28) reads 

_ Jo XlK(i^(^)X2{!;,P) sin^(g_a)sin^(g-/3) 

~ f 27r min(sin^(g-a),sin^(g-/3)) ,^ 

Jo sin^(e-a)sin^(^-/3) "'^ 

= -cos2(a-/3), (31) 

in exact agreement with the quantum mechanical result (see Table I). 
For d = 0, d = 1, and d = 4 we find 

Eia,(3)=-1+'-^^^^^, (32) 

TT 
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1+1 


cos(a- 


-m 


1+1 


sin(a- 




1-1 


cos(a- 


-m 


1-1 


sin(Q;- 





and 

E(a,/3) = _(3-cos^2(a-/?))cos2(a-/3)^ ^3^^ 

respectively. The corresponding results for S'(0) are shown in Fig. 3. 

If ^ 00, @{W - \tn,i - tn,2\) = 1, and P(ri,T2, W) = 1 such that Eq. (28) reduces to^ 

1 r'^'^ 

E{a, /3) = / sign(cos 2(e - a)) sign(cos 2(e - /3))de 

^_^^2|a-/3|mod7r_ ^^^^ 

TT 

Obviously, Eq. (35) does not agree with the quantum theoretical expression E{a^ [5) = 
— cos 2(a — 13). 
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